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The question of whether the Bose-Einstein condensation (BEC) is a necessary condition of the 
Meissner effect is examined. The electromagnetic susceptibility of a charged Bose gas with short- 
range repulsion is studied using the perturbation theory with respect to the repulsive force. With 
decreasing temperature, the Bose-statistical coherence grows, and prior to the BEC phase the sus- 
ceptibility shows a singularity implying the Meissner effect. This means that the BEC is a sufficient, 
but not a necessary, condition of the Meissner effect in the charged Bose gas with short-range re- 
pulsion. 
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I. INTRODUCTION 

The Meissner-Ochsenfeld effect (which from now on 
we call the Meissner effect) is known as a striking phe- 
nomenon showing directly the existence of the Bose- 
Einstein condensate (BEC). Actually, when cooling the 
metal in an applied magnetic field, just at the moment 
when the metal becomes superconducting, the metal be- 
gins to exclude the magnetic field. 

A simple model exhibiting the Meissner effect is the 
charged Bose gas. Before the advent of the BCS model, 
Schafroth gave a clear definition of the Meissner effect us- 
ing the charged Bose gas The BCS model, starting 
from a realistic picture of the electrons in metal, estab- 
lished a realistic model exhibiting the Meissner effect [|[ . 

A common feature of these two models is an appear- 
ance of the off-diagonal long-range order (ODLRO) in 
low temperature such as [J], 

< *(a!)t$(0) >^ f*(x)f(0), |z| — oo. (1) 

To make the problem clearer, a model-independent 
derivation of the Meissner effect from fundamental prin- 
ciples and the ODLRO is desirable. Recently, such a 
derivation was attempted, with the result that the Meiss- 
ner effect is derived from the gauge invariance and the 
ODLRO This result means that the ODLRO is a 

sufficient condition of the Meissner effect. 

For the converse statement that the ODLRO is a neces- 
sary condition of the Meissner effect, however, the prob- 
lem is not so simple. Consider the following counterex- 
ample. In the two-dimensional systems, the order which 
one observes at zero temperature is not the ODLRO, but 
the off-diagonal finite-range order (ODFRO), which falls 
off exponentially at a long distance [3] ■ But its dynam- 
ical response is not an ordinary one. Experimentally, 
the thin-film superconductors exhibit the Meissner ef- 
fect. Theoretically, the two-dimensional charged Bose 
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gas, although it does not show any mathematical sin- 
gularity in the thermodynamic quantities, excludes the 
applied magnetic field at low temperature, which cannot 
be distinguished from the conventional Meissner effect 
in a practical sense |g]. This example poses a question 
of whether the ODLRO is a necessary condition of the 
Meissner effect. In other words, we face another question 
of whether we must have a more general definition of the 
Meissner effect than our present one. 

In this paper, we assert that the ODLRO is not a nec- 
essary condition of the Meissner effect. The Meissner 
effect occurs even at temperatures in which the many- 
body wave function obeying Bose statistics grows to a 
large but not yet macroscopic size (ODFRO) @. 

An objection to this assertion is that one does not 
observe in the superconductors the Meissner effect at 
T > T c . The phenomenon occurring in the superconduct- 
ing metals at the vicinity of T c is, however, not a gradual 
growth of the coherent wave function from a finite to a 
macroscopic scale, but a formation of the Cooper pairs 
from two electrons. Once the Cooper pair forms a com- 
posite boson at low temperature and high density, the 
system does not change through the pre-BEC state, but 
it immediately jumps to the BEC phase Hence, the 
Meissner effect abruptly occurs at T c without any precur- 
sor. In this meaning, the absence of the Meissner effect 
at T > T c in the superconductors does not contradict the 
above assertion. 

In principle, the charged Bose gas is a most elemen- 
tary model of the Meissner effect. But, we do not have 
any experimental example except for the superconduct- 
ing metal. This situation leads us to a preconceived idea 
that the ODLRO is a necessary and sufficient condition 
of the Meissner effect. Although the BEC is a most pe- 
culiar feature of Bose statistics, it does not automati- 
cally mean that all anomalous properties are attributed 
to the BEC. The bosons obey Bose statistics even when 
the Bose condensate is absent. Hence, to what extent 
these anomalous phenomena imply the existence of the 
Bose condensate is, logically, a problem to be consid- 
ered separately. From this standpoint, we can classify 
all phenomena in the dense cold bosons into two classes, 
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the distinguishing factor being the necessity of the Bose 
condensate. The first is a class for which the Bose con- 
densate is a necessary and sufficient condition, and the 
second is a class for which the Bose condensate is only 
a sufficient one, thus being possible even when the Bose 
condensate is absent. 

In real systems, the existence of the interaction be- 
tween particles complicate the above classification. An 
interplay between Bose statistics and the repulsive or 
attractive interaction determines the class of the phe- 
nomenon. In this paper, we study the Meissner effect 
using the charged Bose gas with short-range repulsion 
having the following hamiltonian, 



H 



d A x 



E 2 + B 2 



+ \(d li +ieA li )$(x)\' 



p,p' i 



(2) 



where $ is the spinless boson, and g (> 0) represents 
the short-range repulsive interaction due to the internal 
structure of the bosons (fi — x, y, z). 

As will be explained in Sec.2A, the Meissner effect oc- 
curs when the transverse excitation of the system is sup- 
pressed compared with the longitudinal excitation [111 ]. 
Following Feynman's argument [13], Sec.2B explains that 
this suppression is a fundamental feature of Bose statis- 
tics, being independent of the presence or absence of the 
BEC. In the charged ideal Bose gas (g — 0), where the 
long-range Coulomb force is assumed to be screened by 
the opposite charges in the medium, the suppression of 
the transverse excitation comes only from Bose statis- 
tics, and it exhibits the Meissner effect only in the BEC 
phase. In the charged Bose gas with short-range repul- 
sion (g > 0), however, the collective excitation due to 
the repulsive interaction complicates the situation. As a 
general property of the gas, the longitudinal excitation 
prevails over the transverse one, thus destroying further 
the balance between the transverse and the longitudinal 
excitation. When this effect is added to the suppression 
due to Bose statistics, it is natural to ask whether the 
Meissner effect occurs before the many-body wave func- 
tion grows to a macroscopic scale. 

Sec. 3 studies the Meissner effect by calculating the 
electromagnetic susceptibility of the charged Bose gas 
with short-range repulsion using the perturbation expan- 
sion with respect to the repulsive interaction g. To exam- 
ine a meaning of the BEC for the Meissner effect, we can- 
not assume the Bose condensate from the beginning by 
replacing the operator of zero-momentum particle with a 
c- number [la ]. Rather, it is crucially important to incor- 
porate the Bose-statistical coherence of the many-body 
wave function without assuming the Bose condensate. 
Starting from the normal phase, the perturbation must 
be developed in such a way that as the order of the ex- 
pansion increases, the susceptibility includes a new effect 
due to a larger coherent wave function. Sec. 3 shows an 
occurrence of the Meissner effect prior to the BEC phase 



in cooling the charged Bose gas with short-range repul- 
sion. Sec. 4 discusses the role of the repulsive interaction 
from a different point of view, and compare this Meissner 
effect with the conventional one in the superconductors. 



II. MEISSNER EFFECT IN THE CHARGED 
BOSE GAS 

A. Divergence of the susceptibility 

We apply to the charged Bose gas a weak stationary 
magnetic field derived from a transverse vector poten- 
tial Afj,(q 7 uj)e tuJt . The current density within the linear 
response is, 

a 2 



(J M (g,w)) 



[nm-x (q,^)}Af,(q,uj), (3) 



where n and m is a number density and a mass of the par- 
ticle, respectively. The first term in the right-hand side 
of Eq.(3) is the usual diamagnetic current, and the sec- 
ond one is the paramagnetic term arising from a change 
in the system's wave function induced by A^^q^uj). The 
second term is a transverse part of the Fourier 

component of the current-current response tensor, 

1 f 

Xixu{q,^n) = y / dTexp(iuj n T){T T J fl (q,T)J l/ (q,0)), 

(4) 



where 



j ll ( q , T )=j2[p+$) <f> P <w"^ , (5) 



(h = 1), and defined as, 



Inside of the matter, Eq.(3) is added to the Maxwell 
equation in vacuum, 

-(^-g 2 )^,u,) = ^<4%,a,)>, (7) 

as an additional source term. Hence, one obtains a 
following relationship between A tl (q,u) in matter and 
A e ^(q,Lu) in vacuum, 



An(q,uj) 



Al x (q^) 



4ne ( nm — x (q,&) 



1 - 



U! 2 /C 2 



(8) 



For the static limit, from gauge invariance, one knows 
nm = x L (q> 0) relating the longitudinal current-current 
response tensor y L (q, to) with number density n and mass 
m of the particle [lj]. Hence, one gets, 



^(9,0) 
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(9) 
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Equation. (9) means that, when the balance between 
the transverse and the longitudinal excitation is de- 
stroyed in the static and uniform limit (]im q ->o[x L (q, 0) — 
X T {q, 0)] 7^ 0), the denominator in the right-hand side of 
Eq.(9) diverges as q — > 0, thus vanishing A^(q, 0) and 
H = V x A a,t q — » in matter (Schafroth's criterion of 
the Meissner effect). 

From now on, we express a term proportional to q^q v 
in Xiiv by x^ v such as, 



Hence, Eq.(9) is expressed as, 



x L (q,u) - x T (q,u) 



A^(q,0) = 



A e *(q,0) 



47re 2 /%„(g,0) 



1 



(10) 



(11) 



(1) In the charged ideal Bose gas, Xfif is given by. 

q^q v ^f(e(p))-f(e(p + q)) 
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(b) 



where e(p) is the kinetic energy p 2 /2m of the spinlcss 
boson, and f(e(p)) is the Bose-Einstein distribution. In 
the BEC phase (chemical potential p — 0), /(e(p)) in 
Eq.(12) is a macroscopic number for p — and nearly 
zero for p ^ 0. Thus, in the sum over p in the right-hand 
side of Eq.(12), only two terms corresponding to p = 
and p = —q remain, with a result that, 



q»qu 



(13) 



where n is the number density of the p — Bose par- 
ticles. In view of q~ 2 in Eq.(13), the charged ideal Bose 
gas in the BEC phase exhibits the Meissner effect. In 
the normal phase (/i < 0), however, when the sum over 
p in Eq.(12) is carried out by replacing it with an inte- 
gral, one notices that q~ 2 dependence disappears in the 
result. Hence, in the charged ideal Bose gas, the BEC 
is a necessary and sufficient condition of the Meissner ef- 
fect. In this case, the penetration depth Ao(T), defined in 
Eq.(3) by J^q,0) = - (c/47r)(l/A g(T))A M (g, 0) at q -> 0, 
is A (T) = n (T)- -V m c 2 /47re 2 . 

(2) To analyze the Meissner effect in the charged 
Bose gas with short-range repulsion, instead of 
\im q ^o[x L (q, 0) — X T (<7,0)] ^ 0, we regard a divergence 
of the coefficient of q^q^ in Xiw 

(g,0) at 9 -0 (Eq.(10)), 
as a general definition of the Meissner effect. 



B. Effects of Bose statistics 



The physical explanation of Ym\ q ^\x_ L {q, 0) — 
X T (q, 0)] ^ in the Bose gas dates back to Feynman's 



FIG. 1: A schematic picture of (a) transverse displacement, 
and (b) longitudinal displacement in coordinate space. White 
circles represent the bosons in the initial configuration. 



argument on the scarcity of the excitation in the liquid 
helium 4 [T^. We recapitulate his explanation of how 
Bose statistics affects the many-body wave function in 
configuration space. 

(1) When the gas is in the BEC phase (fi = 0), the wave 
function has the permutation symmetry on a macroscopic 
scale. Consider a collective transverse displacement mov- 
ing particles on a solid straight line in Fig. 1(a) to that on 
a dotted wavy curve (black circles). At first, it seems a 
configuration change on a large scale. This displacement 
is, however, reproduced by a set of slight displacements 
(depicted by short arrows) moving the particle close to 
the wavy curve in the initial configuration to the exact 
positions on the curve. The transverse displacement like 
Fig. 1(a) changes only slightly the local density of par- 
ticles (defined in some volume). Hence, at any part of 
the wavy curve, it is possible to find in the initial con- 
figuration a particle close to the wavy curve. In Bose 
statistics, due to permutation symmetry, one cannot dis- 
tinguish between one particle on the wavy curve which 
comes from the position close to the wavy curve by the 
short arrow and another particle which comes from the 
straight line by the long arrow. Even if the displace- 
ment from the straight to the wavy line is a large dis- 
placement in the classical statistics, it is only a slight 
displacement in Bose statistics (short arrows). Hence, 
the excited state lies in a small distance from the ground 
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state in configuration space. Since the excited state is 
orthogonal to the ground state, the wave function cor- 
responding to the excited state must spatially oscillate. 
Accordingly, the many-body wave function of the trans- 
versely excited state oscillates within a small distance in 
configuration space. Since the kinetic energy of the sys- 
tem is determined by the gradient of the wave function, 
this means that the energy of the transverse excitation 
is not small even at q — 0, leading to the scarcity of the 
low-energy transverse-excitation density. This is a reason 
for x T (g, 0) — > as q — > at low temperature. 

On the other hand, for the longitudinal displacement 
like Fig. 1(b), we see a different situation. Since the longi- 
tudinal displacement changes the local density so much, 
it is not possible to always find in the initial configu- 
ration a particle which is close to a given particle after 
displacement. (Even if it is possible for only a few par- 
ticles after displacement, it is not possible for all parti- 
cles.) Hence, one can not reproduce the configuration 
pattern after displacement only by permutation of the 
initial pattern (a lower pattern in Fig. 1(b) can not be 
reproduced with an upper pattern). This means that the 
longitudinally-displaced configuration is far from the ini- 
tial one in configuration space. Hence, the wave function 
of the longitudinally excited state gradually oscillates in 
a long distance. Applying the definition of the kinetic 
energy to this case, one concludes that the low-energy 
longitudinal excitation is possible. This is a reason for 
X L (l, 0) 7^ at q — ► 0. The above explanation asserts 
that Bose statistics is an essential reason for the Meiss- 
ner effect. 

(2) When the gas is at high temperature (/i <C 0), the 
wave function has the permutation symmetry only in the 
limited area. Accordingly, in the transverse displacement 
like Fig. 1(a), one can not regard the wave function after 
permutation as an equivalent of the initial wave function. 
The displacement which was regarded as a small displace- 
ment in the BEC phase becomes a large one in the normal 
phase. Since the difference between the transverse and 
the longitudinal excitation discussed in (1) vanishes, the 
Meissner effect disappears at high temperature. 

(3) From our viewpoint, a remarkable state of the Bose 
gas lies at the vicinity of the BEC transition temperature 
in the normal phase (fi < 0) [l5l |. The coherent many- 
body wave function grows to a large size but not yet to a 
macroscopic one. In such a situation, whereas the mech- 
anism in (2) works for the large displacement extending 
over two different wave functions, the mechanism in (1) 
works for the small displacement within a single wave 
function. 

As shown in Eq.(13), the BEC is the necessary and 
sufficient condition of the Meissner effect for the charged 
ideal Bose gas. This conclusion is, however, applicable 
only to a simplified model like the charged ideal Bose 
gas. An inclusion of the repulsive interaction will change 
the situation, which will be studied in Sec. 3. 
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FIG. 2: The first-order Feynman diagram of the current- 
current response tensor among the one-particle excitations 
of the repulsive Bose gas in Eq.(14). (The solid and the dot- 
ted lines represent the boson and the repulsive interaction, 
respectively.) Particle exchange between the tensor and the 
medium by Bose symmetry is explained in four steps from (a) 
to (d). 



III. FORMALISM 

To examine the normal Bose gas at the vicinity of the 
BEC transition temperature, we can not use the con- 
ventional method of assuming the Bose condensate by 
replacing the operator of zero momentum particle with 
a c-number from the beginning. Consider the current- 
current response tensor of the charged Bose gas at the 
vicinity of the BEC transition temperature as illustrated 
in Fig. 2. The Xftv m coordinate space includes the fol- 
lowing form, 



(G\T t J^(x,t)J v (0,0)\G) 



(0\T T J IM (x,T)J v (0 ) 0)exp 


\~f 

Jo 


dr Hj (t ) 


10) 


(0\exp 


-/ drflj(r) 
Jo 


10) 



'(14) 



where \G) is a ground state of J2 p € (p)^I^p + 

Ep, P ' E 9 ®i- q ®l>+ q ®p'®P- The excitation energy due 
to the long-range Coulomb force is too high to affect the 
dynamics of the bosons in the Meissner effect [16j. Hence, 
the photon-boson vertex appears only at an initial and a 
final vertex of the x^v-, and Hi(t) in Eq.(14) represents 

9E P , P >EAl- q ®l> +q %'% ( dotted lines in Fi s- 2 )- In 

the ground state of Bose gas, due to the repulsive inter- 
action g contained in exp(—i J Hi(t)<1t) of Eq.(14), one- 
particle excitations frequently occur as illustrated by an 
upper bubble with a dotted line in Fig. 2(a). The current- 
current response tensor of such a medium is depicted by 
a lower bubble in Fig. 2(a). (The black and the white 
small circle in Fig. 2 represents the vector and the scalar 
vertex respectively.) 
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When \G > is the ground state at the vicinity of the 
BEC transition temperature in the normal phase, the 
perturbation must be developed in such a way that as 
the order of the expansion increases, the susceptibility 
gradually includes a new effect due to Bose statistics. An 
important feature of the response to the photon by the 
large coherent many-body wave function is that bosons 
appearing in the response tensor also participate in the 
coherent wave function of the medium. (The response 
tensor and the medium form a coherent wave function as 
a whole.) Hence, we must seriously consider the influence 
of Bose statistics on the graph like Fig. 2(a). When one 
of the two bosons in the lower bubble and one of the 
two bosons in the upper bubble have a same momentum 
(p = an d another bosons in the lower and the upper 
bubble have another same momentum (p + q — p' + q'), 
a graph made by exchanging these two particles must 
be included in the expansion of Xnv Such a graph is 
obtained in two steps as illustrated in Fig. 2. When two 
boson lines with p and p'{= p) is exchanged in Fig. 2(a), 
a square appears in Fig. 2(b). Further, when two boson 
lines withp+g andp+<j'(= p+q) is exchanged in Fig. 2(c), 
a graph linked by the repulsive interaction appears in 
Fig. 2(d). The resulting current-current correlation tensor 



has a form elongated by the density-density correlation. 
A contribution of Fig. 2(d) to x^v is given by, 

x r_ Mv))- m± m\ (15) 

u) + e(p) - e(p + q) 

With decreasing temperature, the coherent wave function 
grows to a large size, and the particle exchange due to 
Bose statistics like Fig. 2 occurs many times. Hence, one 
can not ignore the higher-order term x$ corresponding 
to the larger coherent wave function. Particularly im- 
portant is the coherent wave function consisting of p = 
bosons. Correspondingly, a quantity of most interest 
which contribute to x^v in Eq.(lO) is a term proportional 

to q^q v in xffiili 0), the p = component of which has 
the following form; 

X^(^)= qj f L 9 n F P {q) n+ \ (16) 

where 



F p {q) = exp(/3/x) 



(1 - exp^/i))- 1 - (cxp(/?e( g )) - exp(/^)V_ 



(17) 



r 



Fp(q) is a positive monotonously decreasing function of 
q 2 which approaches zero as q 2 — > oo. As fi — ► 0, l-F^t?)! 
increases at any q. 

The general definition at the end of Sec.2A says that, 
when a power series in g 



n+l 



2 2 

n=0 n=0 

(18) 

diverges at q — ► 0, the Meissner effect occurs. At high 
temperature in normal phase <C 0), a small F@(q) 
guarantees the convergence of X/^(<7>0) (Eq.(18)). With 
decreasing temperature (fi — > 0), however, a gradual in- 
crease of Fp(q) makes the higher-order term significant, 
finally leading to the divergence of x^(q^)- A conver- 
gence radius r c of the power series X^^Li a n,x n is given 

by 



= lim M 1 /", 



(19) 



(Cauchy-Hadamard's theorem). Applying this theorem 
to Eq.(18) (a n = F fj {q) n+1 ), one obtains l/r c = F p {q). 
Hence, as a convergence condition: < g < r c , we have 



\gFp(q)\ < i- 



(20) 



Generally, for the Meissner effect to occur, the balance 
between the longitudinal and the transverse excitation 
must be destroyed on a macroscopic scale (q — > 0). To- 
gether with the Bose-statistical coherence, the repulsive 
interaction enhances this tendency, finally violating the 
condition of Eq.(20). 

An expansion form of Fp(q) around q 2 = is given by 



Fp{q) 



/3exp(/3^) 
(1 - exp(/3 M ))2 



p 2 (1 + expjfiii)) 

- T exp W (i _ cxpW)3 £(g ) + . 



(21) 



At high temperature, Eq.(20) is satisfied, and Xnv{q,fy 
has the following form, 



X^(9,0) 



2 i-gFp{ q y 



(22) 



In cooling the gas, |-F/3( g )| increases monotonously, and 
the convergence condition is first violated at q — when 
|£i.F>(0)| = 1, that is, using Eq.(21), 



5 /3 = 4sinh^ ( ^ 



(23) 
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From now, we call To an onset temperature of the Meiss- 
ncr effect satisfying Eq.(23). In view of Eq.(23), wc con- 
clude that the charged Bose gas with short-range repul- 
sion (g > 0) shows the Meissner effect prior to the BEC 
phase (fj, < 0), thus leading to T BEC < T . At T = T , 
substituting Eq.(21) in Eq.(22) with Eq.(23), we get 



m 



sinh q 2 



(24) 



which shows the Meissner effect in Eq.(ll). 

Compared with Eq.(13), the number density no of 
p = boson which is a macroscopic number at Tbec 
is replaced by 0.5/ sinh |/3/x| in Eq.(24). Accordingly, 
the penetration depth has a form such as A(T) = 
\J sinh \j3jj\y/ mc 2 / Aire 2 . Since /x < at T , A (To) is larger 
than Xo(Tbec) m the ideal Bose gas, because at To only 
p = boson whose number is large, but not yet macro- 
scopic is responsible for this Meissner effect. Hence, at 
the onset temperature To, the magnetic field penetrates 
deeply into the gas. (It is in contrast with the abrupt ex- 
clusion of the magnetic field, that is, a very small A(T C ) 
in the superconductor.) 

(a) In the weak-coupling case (g ~ 0), we can make 
a rough estimation of the condition of Eq.(23). Using 
iio(T) of the ideal Bose gas for it(T) in Eq.(23) such as 



/?mt) 



_ (2M\ 2 T BEC 



\ Tbec J 



1.5 



(25) 

we can determine T as a solution of simultaneous equa- 
tions (23) and (25). Figure. 3 shows a phase diagram in 
which the onset temperature T of the Meissner effect is 
plotted at a given strength of the repulsive interaction g. 
(Since only the relative strength of g to the particle mass 
has a meaning, g/ksTsEC is plotted.) Figure. 3 confirms 
that for the charged ideal Bose gas (g = 0), the Meiss- 
ner effect occurs just at T = Tbec, but it shows that 
the charged Bose gas with short-range repulsion (g ^ 0) 
begins to exclude the applied magnetic field prior to the 
BEC phase in cooling. (Although the repulsive interac- 
tion suppresses the Tbec, it does not change the overall 
feature of Fig. 3.) The BEC is not a necessary condi- 
tion of the Meissner effect for the charged Bose gas with 
short-range repulsion. 

(b) At high temperature (T/Tbec I), the repul- 
sive interaction which is necessary for the Meissner effect 
to occur in Eq.(23) increases as Texp[a(T/Te£;c) 2 ]- In 
view of the exponential factor, the Meissner effect does 
not occur at high temperature in real materials, thus not 
affecting the essential feature of the Meissner effect, a 
phenomenon in low temperature. 

(c) When cooling the system further below To in Fig. 3 
at a given g , the x^iq, 0) (Eq.(18)) diverges not only at 
q = 0, but also for a region < q < q c {T), where q c {T) is 
determined by a condition \gFp(q c (T))\ = 1 at a given T. 
For <7 C (T) < q, the x^(q, 0) converges even below T . To 
obtain the overall feature of x MI y(g, 0) below T , we get 




FIG. 3: A phase diagram showing the onset temperature of 
the Meissner effect To at a given strength of the repulsive 
interaction g. Both quantities are normalized by Tbec and 
ksTBEC, respectively. 



a concrete form of Xv-v by expanding 1 — gFp{q) in the 
right-hand side of Eq.(22) with respect to q 2 - q c (T) 2 . 
Hence one gets the following result: 



XV(<7,0) 
( oo, 



< 



25 



-q^q v 



dlnFp(q) 
8q 2 



(q 2 -q c (Tf) 



0<q< q c (T), 
q c (T) < q. (26) 



Equation. (26) in Eq.(ll) means that the charged Bose 
gas with short-range repulsion excludes not only the uni- 
form applied magnetic field (q = 0), but also the spatially 
modulated magnetic field H(q) with < q < q c (T). 

With decreasing temperature, the number of particles 
with p ^ which satisfies the divergence condition of 
X/if(9)0) gradually increases. Hence such particles par- 
ticipate in excluding the applied magnetic field, leading 
to a gradual decrease of the penetration depth. 

The q c (T) reminds us of the coherence length £ in the 
superconductor, because q c {T) and have a similar 
property: As they become larger, the Bose-statistical co- 
herence grows. (For more comparison, see Sec. 4.) Here, 
we discuss the temperature dependence of q c (T). 

(1) Just below T), q c (T) in Eq.(26) is small, and we can 
obtain the temperature-dependence of q c (T) as follows. 
We expand Fp(q) in the condition of \gFp(q)\ = 1 with a 
small q c (T) as in Eq.(2I). Taking only the first and the 
second term in the right-hand side of Eq.(21), we get an 
approximate form of q c (T) as 



q c (T) 2 = 4mfc B Ttanh 



\MT)\ 



1 



4fc B T 



■ sinh 



2 
(27) 

At T = T , q c {T) is zero because of Eq.(23). At T < T , 
q c {T) is expanded with respect to Tq — T. In the case 
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of weak repulsive interaction, using Ho(T) (Eq.(25)) as 
fi(T), one gets 



q c (T) = 2\j Amk B T tanh 
where 



\PMT)\ 



(28) 



A = 1 + 1.08(a 2 - 0.5a 05 - 0.5a~ 1 ) ] j^^+l, 



(29) 



and a = T /T B ec- 

(2) When the decreasing temperature reaches T = 
Tbec, Fp{q) (Eq.(17)) diverges at all q, thus leading to 
a divergence of Xnv{<li 0) at all Q- Hence q c (T) in Eq.(26) 
goes to oo in the BEC phase. 



IV. DISCUSSION 

(1) In this paper, we explained the mechanism in 
which the repulsive interaction destroys the balance be- 
tween the longitudinal and the transverse excitation on a 
macroscopic scale. This explanation can be viewed from 
a different point of view. 

Generally, the influence of the repulsive interaction on 
the dynamical properties of the boson systems differs 
from that on the thermal properties. For the thermal 
properties, the repulsive interaction generally suppresses 
the growth of the Bose-statistical coherence. For exam- 
ple, it is known that due to the strong short-range repul- 
sion, not all helium 4 atoms participate in the BEC even 
at zero temperature (depletion effect). 

For the dynamical properties such as the superfluid- 
ity, however, the repulsive interaction generally enhances 
singular properties due to Bose statistics. For the ideal 
Bose gas, the particle with non-zero momentum behaves 
like a free particle (except for the statistical constraint) 
to the external dynamical perturbation. For the repul- 
sive Bose gas, however, the particles are likely to spread 
uniformly in coordinate space due to the repulsive force. 
This feature makes the particles with p^O behave sim- 
ilarly with other particles, especially with the particle 
having zero momentum. If it behaves differently from 
others, a resulting locally high density of particle raises 
the interaction energy. This is a reason why all particles 
seem to show singular dynamical behavior even at tem- 
perature in which not all particles participate in the BEC. 
(This fact corresponds to the well-known property that 
the one-particle energy spectrum changes from p 2 /2m to 
sp due to the repulsive interaction, enhancing the stabil- 
ity of the superfluid flow to the thermal disturbance by 
raising the energy barrier around p = 0.) 

The Meissner effect is a dynamical response of the Bose 
gas to the external perturbation (the applied static mag- 
netic field) as well. The boson excited from the p = 
to the p ^ state does not respond as a free particle 
to the magnetic field, but it excludes the magnetic field 



cooperatively with the p = bosons by the repulsive 
interaction. Hence, even when the Bose-statistical co- 
herence grows to a large but not to a macroscopic scale 
(ODFRO), the charged Bose gas with short-range repul- 
sion shows the Meissner effect. In this meaning, as well 
as the superfluidity, the Meissner effect is another exam- 
ple of the mechanism in which the repulsive interaction 
enhances the singular dynamical response of the boson 
system. 

(2) As discussed in Sec,l, the Meissner effect occurring 
in the charged Bose gas with short-range repulsion dif- 
fers from that in the superconductors in some respects. 
Comparing two different Meissner effects is useful for a 
deeper understanding of them. 

Let us replace the boson operator <E> p in the right-hand 
side of Eq.(5) with the fermion operator ^ p ,cr, and use 
such a Jft(q, t) in Eq.(4). By factorizing the resulting 
{T T J l _ L {q 1 r)J v {q, 0)) into a product of the two anomalous 
fermion Green's functions, we get 

i n \ 1 f d 3 p . q. , q. 
X^(<?>°)=2 J ^f(P+2'^ p+ 2' l/ 

(30) 

J-'pJ-'p+q J -Up + ^p+q 

where e p = p 2 /2m — \i, A is the order parameter, and 



E p = ^Je 2 + A 2 . Performing the integral over p in 
Eq.(30) yields 



X T (q,Q) = (nm-a\A\ 2 )-bq 2 



(31) 



where nm = x L (q> 0)j and a and & is a parameter inde- 
pendent of q 2 . Substituting Eq.(31) in the definition of 
Eq.(10), one obtains 

X IM ,{q,0)=q li q v (^- + bj, (32) 

which means that only below T c (|A| ^ 0), X^g, 0) 
shows the mathematical singularity like q~ 2 , leading to 
the Meissner effect. 

As well as the charged ideal Bose gas, the BCS model 
does not show the mathematical singularity in X/^feO) 
prior to the BEC, but after the BEC, only the first-order 
term in the perturbation expansion of Xnv(li 0) is enough 
to show the q~ 2 singularity. In the BCS model, the par- 
ticle interaction is assumed to be completely diagonal- 
ized by the formation of the Cooper pair, thus leaving no 
residual interaction between the pairs. This fact makes 
the next-order calculation of the susceptibility unneces- 
sary. On the other hand, the charged Bose gas with short- 
range repulsion shows another type of singularity prior to 
the BEC in the normal phase, if the perturbation expan- 
sion of Xhv{q,Q) is summed to infinite order. 

In the superconductors, two different kinds of the co- 
herence length are known. The first is £o appearing in 
the Pippard nonlocal electrodynamics, and the second is 
£(T) appearing in the Ginzburg-Landau equation. They 
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are not the same quantities, but both quantities repre- 
sent the smallest size of the wave packets which the su- 
percurrent carriers can form, thus reflecting an extent 
of transformation of individual particles from fermion to 
boson. Due to the composite-boson nature underlying 
the Cooper pairs, the type-2 superconductors show the 
complicated nonuniform response to the uniform applied 
magnetic field. 

On the other hand, the response of the charged Bose 
gas with short-range repulsion to the magnetic field is dif- 
ferent, because it consists of the elementary boson. (£ 
or £(T) corresponds to the particle radius in this case.) 



The critical wave number q c (T) defined as Eq.(28) is an 
inverse of another type of the coherence length. The 
l/q c (T) represents a characteristic length of the spatial 
variation of magnetic field. Generally, compared with the 
uniform field case, screening of the spatially modulated 
applied magnetic field needs the stronger Bose-statistical 
coherence in the gas. When the applied field varies with 
a length larger than the value of l/q c (T), it is screened by 
the charged Bose gas with short-range repulsion. The in- 
crease of q c {T) in Eq.(28) means the growth of the Bose- 
statistical coherence in cooling the system. 



[1] M.R.Schafroth, Helv.Phys.Acta,24, 645 (1951). 
[2] M.R.Schafroth, Phy.Rev, 100, 463 (1955) 
[3] J.Bardeen, L.N.Cooper and J.R.Schriefer, Phy.Rev, 108, 
1175 (1957) 

[4] O.Penrose, Philos.Mag, 42, 1373(1951), O.Penrose 
and L.Onsager, Phy.Rev, 104, 576 (1956), C.N.Yang, 
Rev.Mod.Phys, 34, 694(1962). 

[5] G.L.Sewell, J.Stat.Phys,61, 415(1990) 

[6] H.T.Nieh, Gang.Su and Bao-Heng Zhao, Phys.Rev.B, 51, 
3760(1995) 

[7] T.M.Rice, Phys.Rev.A, 1889(1965), P.C.Hohenberg, 

Phys.Rev, 153,383(1967) 
[8] R.M.May, Phys.Rev, 115,254(1959) 

[9] Note that ODFRO is a different concept from the diag- 
onal long-range order (DLRO) found in various ordered 
state. 

[10] The fluctuation occurring at the vicinity of T c in the su- 
perconductors is smaller than that appearing in the liquid 
helium 4 close to Tbec- 



[11] As a review, see G.Baym, in Mathematical methods in 
Solid State and Supcrhuid Theory (ed by R.C.Clark and 
G.H. Derrick), 121 (Oliver and Boyd, Edingburgh, 1969) 

[12] R.P.Feynman, Progress in Low Temp Phys. 1, (cd 
C.J.Gorter), 17 (Amsterdam, North-Holland 1955) 

[13] It is well known that the charged Bogoliubov model, in 
which 3>p=o is replaced by the c-number in the hamil- 
tonian of Eq.(2), shows the Meissner effect like the su- 
perconductors (Eq.(13)) at the first order of g. A theory 
without such a replacement is a purpose of this paper. 

[14] See, for example, P.Noziere and D. Pines, The Theory of 
Quantum Liquids 1, (Benjamin, New York, 1965) 

[15] S.Koh, Phys.Rev.B, 64, 134529(2001) 

[16] The energy spectrum of the charged Bose gas with long- 
range Coulomb repulsion has been studied by many au- 
thors. S.K.Ma and C.W.Woo, Phys.Rev, 159, 165(1967), 
and references therein. 



